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ABSTRACT 

Time domain dipole antenna array studies conducted at the Electromag- 
netic Effects Laboratory* of the Harry Diamond Laboratories for application 
to large scale nuclear electromagnetic pulse (EMP) simulation are described. 
Analytically '.nd experimentally derived EMP environments for a single 
dipole antenna, are compared with corresponding data for multiple element 
broadside and collinear dipole antenna arrays. The multiple element concept 
is studied as a technique for eliminating the need for extremely high volumes 
in the single element type nuclear EMP simulators by substituting the super- 
imposed electromagnetic fields radiated by several lower voltage sources to 
obtain the - bivalent EMP environment. Experimentally measured input 
impedances «_re compared with the results from classical antenna theory to 
determine the effects of mutual coupling between elements in the array 
configurations. The aperture size of the array is shown to be an important 
factor in determining the pulse shape characteristics of the EMP environment 
in the near and far fields. Results of the studies define specific advantages 
and limitations in the dipole antenna array techniques for nuclear EMP 
simulation. 


* Formerly with USA Mobility Equipment R&D Center, Ft. Beivoir, Vo. 22060 
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1. INTRODUCTION 


The RES (Radiating Electromagnetic Simulator) type antennas have been used success- 
fully as nuclear electromagnetic pulse (EMP^ simulators, ^hese simulators are dipole an- 
tonnas that are driven in the time domain by a MARX generator located at a blconical center 
food soction. A radiating dipole antenna is formed by extending cylindrical arms froih the 
biconos. Tho major limitation of these types of simulators has been the relatively loW volt- 
ages achieved to date for a single source pulser. The voltages have generally been less than 
2 MV for operational dipole systems. A source in this range cannot meet the large area 
covorago requirements such as those that exist for the SAFEGUARD site. It is therefore 
necessary ,either | to increase the voltage for a single source or to superimpose the electro- 
magnetic environments of several sources to provide the required threat-level coverage. 

Much work has been done in antenna theory on continuous wave, (CW) dipole array sys- 
tems. Such arrays have proved their usefulness in many applications. It is of interest, there- 
fore, to examine the properties of these arrays for application to time-domain EMP simula- 
tion. In theory, one can go from a frequency domain to time domain by the Laplace or Fourier 
transformations. These are linear operations through which some of the useful properties of 
CW arrays can be expected to carry over to the time domain. . 

The work presented in this report was conducted at the Electromagnetic Effects Labora- 
tory HDL. Whenever possible, theoretical results have been verified by experiment. Several 
time domain dipole arrays were studied in order to determine the possibility 'of application of 
such arrays to high field level nuclear EMP simulation. 


i 

2. THE SINGLE-DIPOLE ANTENNA i 

» 

Consider first the single horizontal dipole antenna to establish a correlation between 
theoretical predictions of the radiated EMP environments and the experimental results ob- 
tained irom the models studied in this report. The early time EMP fields of a single-dipole 
antenna can be predicted with the Hertzian dipole equations which, in the frequency domain, 
are given below. 1 ; 


(to) = A (a.) e"! 1 * 




a> 


- j k r f 2 Z 0 2 


E r (w) = A (u) c‘» kf 1+ = — cos P 

V r 2 j we r 3 / 


( 2 ) 


E 0 (w) = A (&•) e 


•)kr [ j W M + + 


s in 6 i 


( 3 ) 


r 3 j (. £ r 3 


where \{u) * I(w)h/4rr, I (a.) is the antenna current, h = equivalent length of the antenna, 

z„ - k - w/m « . w ■ rad/aec, t i and e are the permoublllty and permittivity of free space, 


^Romo, S., Whinntry, J. R., Fltld* ond W ov«» In Modom Radis. Soeond Edition, John Wi toy & Sons, Inc., N.Y. 
1964, p. 498. 
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respectively. The geometry for these and subsequently discussed components is shown In 
figure 1, The time domain expression for the antenna current I(t) is represented by: 


I(t) = K \(1 + a, e"* 1 * - (1 + a 2 t) e } 


where K = h/4ir determines the amplitude of the current, and and a 2 , respectively, de- 
termine the faU and rise times of the current. A plot of I(t) is shown in figure 2. The value 
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Figure 2. Normalized input current pulse shope 
computed from equation (4). 


for the factor K is determined experimentally from either a direct current measurement on 
the antenna, or from a measurement of in the radiation field of the antenna. 


2.1 Time-Domain Expressions for the Electromagnetic Fields 

A time-domain solution for equations (1), (2), and (3) is obtained by the following procedure. 
The j« variable is replaced by the Laplace variable s = j w and the Laplace transformation is 
carried out. The results are convolved with equation (4) to obtain the time-domain expressions 
for the dipole fields: 



where K = h/4n is assumed independent of frequency. This assumption requires that the di- 
pole impedance remains constant over the pulse length of interest, and limits the solution to 
the early time regime. That is, a time < (2h 1 /c) in figure 1, where the dipole antenna does 
not receive ground reflections which affect impedance of the antenna. In all cases under 
study, the antennas are located at a height h ; above the ground so that the rise time is less 
than the ground interaction time. We therefore expect a good representation of the peak fields 
from equations (5), (6), ard (7). 
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Ground reflections iu the radiation field alter the free-spare waveshape and must be taken 
into consideration when comparing theoretical and measured results. In the high-frequency 
(early time) range, the ground acts as a good conductor and thus the theory of images can be 
applied by considering the ground as a perfectly conducting plane. The reflected (image) fields 
in the presence of the perfectly conducting ground plane have the same form as equations (5), 

(6), and (7) but arrive at the point Q with a time delay corresponding to the path length differ- 
ence transit time p/c between the incident and image fields, where p is the path length differ- 
ence and c is the velocity of propagation. Orientation ox the incident and image fields in figure 
1 assure that the boundary conditions are satisfied at the conductor surface (h 2 = 0). We also 
note that the Fresnel reflection coefficient assumes the correct limiting values of +1 and -1 for 
vertical (electric field parallel to plane of incidence, £ x ,E y ) and horizontal (electrical field per- 
pendicular to plane of incidence, E t ) polarizations, respectively. 

The resultant fields in the presence of the ground plane are obtained by transforming 
equations (6), (7), and (8) from spherical to rectangular coordinates, and by combining the 
incident and image components with proper time delay and orientation. This procedure gives 
the x, y, z components of the EMP environment, which corresponds to the measured fields. 

The equations are given below, where the "prime" symbol denotes the imt-ge fields, and the 
subscript T represents the total, or resultant, field intensity. 


H r -(-H 4> sina-H^Fin<K) J - (H^ cos a - cos f) J l 8 ) 

Ext = cos a ~ E(J cos >K) cos n - CE r cos a - E ' cos *P) sin 0 ^ 

E yX = (E r sin a + E' sin V) sin J - (E^ sin a * Ej sin V) cos 

E,t = (E 9 - E^) sin - 1 , (E, - E;> cos <? < n > 


2.2 Numerical and Experimental Result s 

Equations (8) through (11) have been programmed on a digital computer and numerical 
results have been obtained for comparison with the experimental results. The factor K in 
equation (4) has been determined experimentally from the radiated H,. to be 38.42 amperes 
(A). Figure 3 shows the experimental configuration of the single-dipole antenna. A 30-kV 
power supply was used as the source to charge a pulser, which in turn fed a 30-ft cable sec- 
tion of 100-n impedance to the center section of the antenna. Three orthogonal components of 
the electric and three orthogonal components oi the magnetic field were jnapped at several 
ranges on and off the centerline (z = 0 plane) of the mapping area. Representative results are 
shown in figures 4 and 5, comparing the theoretically and experimentally derived field com- 
ponents. The curve marked with the z symbol in the magnetic field plots represents the free 
space incident field The other plots correspond to the measured data shown in the oscil- 
loscope pictures. The peak-measured field strengths are designated by each oscilloscope 
picture throughout this report. 
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Figure 3- The single-dipole antenna experimental configuration. 


2.3 Evaluation of Results 


r 
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The components shown in figure 4 correspond to the TEM components cf a plane wave. 
These are designated PRINCIPAL COMPONENTS. Off the centerline, the environment be- 
comes nonplanar with the addition of the E xT and the E yX components. These components 
are not relatable to a horizontally polarized threat EMP and will be designated NON- 
PRINCIPAL COMPONENTS. 

To correlate the simulator EMP environment with a particular threat EMP environment 
of interest, it is necessary to know the equivalent free field EMP environment radiated by the 
simulator. The EMP threat wave-shapes and amplitudes can then be directly compared with 
the simulator output. In general, the field-measurement sensors and the simulator itself are 
operated in the proximity of the ground. Because of this, the free-space environment cannot 
be measured directly. The theory of plane wave reflection at a smooth ground surface (or a 
perfectly conducting plane) can be used to obtain the free field from the measured tangential 
magnetic field (l.e., the magnetic-field component parallel to the ground surface). The tan- 
gential component of the incident magnetic field in figure 1 is the ’lx" component of H^. Boun- 
dary conditions require that the tangential component of (the reflected field) be identical 
in sign to the tangential component of the incident field so that the tangential magnetic field 
will net vanish at the reflecting surface. At points above the surface, the two fields add 


A 
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constructively (since both are of the same polarity) with a time shift equal to the path length 
difference transit time p/e between the incident and reflected fields. If the point of field 
measurement is close to the ground surface, p/c is less than the rise time of the free space 
incident pulse, and the resultant field H tT preserves the waveshape of the incident free field. 

If such is the case, the incident and reflected tangential components add with very little time 
delay, and the only noticable, difference in waveshape between the resultant field, H jiT , and the 
free field, H^ ? is a slight bump on the rising portion pf H iT * This can be verified in figure 4 
by comparing the measured and computed waveforms of H, T with We also note that the 
magnitudes of the computed and measured values of H agree closely, indicating that the as- 
sumption of perfect ground conductivity is justified in tne early-time regime fc v the compo- 
nents that interfere constructively with ground reflections. The components that interfere 
destructively with the incident fields upon reflection, are more sensitive to the infinite ground 
conductivity assumption, as evidenced by the smaller computed values of H yT and E com- 
pared with the measured results. These components vanish at the surface of a perfect con- 
ductor, but are greater than zero at the surface of an imperfect conductor, such as the earth. 
Hence, the measured values .can be expected to be greater than the computed values near the 
surface of the ground. The same trends appear for measurements of H yr and E JT at a posi- 
tion off the centerline as shown in figure 5. The measured waveshape of H , however, differs 
significantly off the centerline from the computed pulse shape, indicating a discrepancy between 
theory and experiment. The difference is due to radiation off the unterminated ends of the an- 
tenna used in the experiment, which is not included in the theoretical consideration. The adr 
ditional radiation results from an abrupt current reversal at the end .of each antenna arm, 
causing a magnetic field of opposite polarity to radiate which interacts with the principal mag- 
netic 'field and changes the measured waveshape. The nonprincipal electric-field components 
agree well in pulse shape and amplitude. Therefore, for positions off the centerline, the theory 
is useful in predicting the relative amplitude of these components compared with the principal 
fields. The equations presented in this section have provided; a method for predicting the EMP 
environment radiated by a pulsed-dipole antenna in the presence of the ground surface. The 
computations have been supported by experimental results. Hie subsequent sections of this 
report extend the methods- developed for the single-dipole antenna to arrays of such antennas 
to determine their possible value as nuclear EMP simulators. 1 ! 


3. FIVE-ELEMENT BROADSIDE ARRAY 

, The field strengths produced by a single-dipole antenna nuclear EMP simulator are limited 
by the magnitude of the current and voltage that can be supplied to the antenna by a single- 
pulse source. This section considers a technique of increasing the simulated EMP environ- 
ment by superimposing the fields of several simultaneously activated antennas which are dis- 
tributed in an array to form a radiating aperture. The five-element broadside dipole array 
model was built to .study this concept. Illustrations of tbe array are given in figures 6 and 7. 
The antennas are driven in parallel from the 30-kV source by five 30-ft-long 100- n cable con- 
nections. 1 This configuration provides, each antenna with a current equivalent to the single- 
dipole current discussed in section 2. As a result, the array aperture contains five times the 
amount of current compared to that of the single source. 


3.1 Superposition of Electromagnetic Fields 

, The principle of superposition will be applied to compute the EMP environment radiated 
by the array. First, it is pointed out that although five times the single- source current exists 
in the array, the radiated environment is not necessarily assured to be five times as large. 
This is a consequence of two important factors. 
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Figure 6. Five-element broadside dipole array geometry. 
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(1) The array is a distributed source , extending over several pulse rlsetimes in its di- 
mensions. The fields computed at some distant observation point result from a retarded time 
superposition of the individual antenna environments. The retarded time refers to the differ- 
ent arrival time for the EMP wavefront from each source to the point of observation due to 
different path lengths between each antenna and the observation point. Thus, the addition of 
the fields from each source is spread out over a time equal to the maximum transit time dif- 
ference between two extreme antennas of the array and the observation point. This results in 
a field magnitude smaller than that produced by a set of simultaneously superimposed pulses. 
Consequently, fields from the five-element array are expected to be less than five times that 
of a single source at finite ranges of observation. At an infinite distance of observation from 
the array, any finite aperture will appear as a point source and the retarded time effect will 
not be observed. However, for EMP applications, energy limitations preclude operation in 
the extreme far field of the array. 
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(2) The second important factor to consider is the possibility of losses occurring from 
mutual coupling between antenna elements. The radiating electric and magnetic fields of each 
antenna eventually strike every other antenna in the array, and some type of interaction (mu- 
tual coupling) must occur. From elementary electromagnetic theory, one expects the coupling 
to be in the form of induced currents which are smaller than, and tend to oppose, the initial 
source currents on the antennas. Thus, in addition to the retarded time effect, mutual coup- 
ling can be expected to further reduce the magnitude of the EMP environment radiated by an 
experimental array. As a starting point only the effects of retarded time superposition will 
be included in the array field equations, and the study of mutual coupling effects is considered 
in a later section of this report. 


3.2 Time-Domain Expressions for the Electromagnetic Fields 

Ideally, each antenna radiates the electromagnetic fie ,ds described by equations (5), (6), 
and (7). The equations corresponding to equations (8) through (11) for the single antenna are 
expressed for the array in terms of the retarded time summation over the number of antennas, 
N = 5, in the array: 


N 

H *T = -2] (H ^ Sina '' +H ^ SinS ’n) < 12 > 

n-1 

N 

Hyr = (H * n COS a " ‘ H *" Cos S '" ) < 13 ) 

n*l 


E *T r T {[E 0r> COS a n - E *n CO ' J V COS " ( E rn COS a n “ E ™ COS * J S » n < 14 > 


n*l 


E yT = {(E rn Sin a n + E ™ sin V sin 9 n ~ S ‘ n + E ffr. sin COS <9 „> < 15 > 

n*J 


N 

V -1 


E «T =2_, {(E ^n - E 9n) Sin * C E rn “ COS 


(16) 


where the "prime" symbol denotes the reflected (or image) field in each equation. The n sub- 
scripted field components contribute to the summations for the resultant fields only after the 
propagation time from the nth source of the point of observation Q has elapsed. Figure 1 shows 
the angles and field components for a typical element of the array. Equations (12) through (16) 
apply for any cumber N of mutually parallel dipole antennas which are located at arbitrary 
positions parallel to tho ground surface. The nth field component, its corresponding retarded 
time, and the angle a n , V n , and 9 r are computed by using the geometry illustrated for a typi- 
cal source and observation point in figure 1. The equations have been programmed on a digital 
computer to provide numerical results. 
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3.3 Numerical apd Experimental Results 


Field measurements were conducted on the experimental model for comparison with the 
theory. A typical set of results is shown in figures 8 and 9 for the observation point at 40 ft 
on and off the centerline in the test area. As was the case for the single-dipole antenna, the 
computed plot represents the free space radiated magnetic field of the array. The T sub- 
scripted components correspond to the ground interacted fields (i.e., the incident free spate 
plus the giuund-reflected environments) which represent the measured fields shown in the 
oscilloscope pictures in figure 8. Again one observes that the measured tangential magnetic 
field H xT compares closely in waveshape with thetheoretical H xT and plots, and repre- 
sents a reliable means of determining experimentally tlie free-space output of the source. 

In comparing magnitudes, the theoretical value for H XT is greater than the measured value. 
This indicates that the omission of mutual coupling from the computations has affected the 
results. The reflection at 35 nsec corresponds to the mismatch at the unterminated ends of 
the antennas. 

Figure 9 shows the measured environments off the centerline of the test area. In this 
region, the nonprincipal components are significant in magnitude, just as they were in the 
single antenna case, and the environment is nonplanar. 


3.4 Evaluation of Results 


A study of the computed electric free field indicates that the peak free field for the five- 
element array is approximately five times greater than the free field of the single antenna. 
Retarded time superposition has not affected the magnitudes appreciably for the particular 
array dimensions under study in this section. This becomes apparent when the path length 
difference transit time is computed between the extreme elements I and V In figure 6. The 
time difference is 1.2 nsec which is less than the 4 nsec 10 to 90% free field pulse rise time 
of each antenna. Consequently, the retarded time degradation occurs within a time span of 
1.2 nsec and is beyond the resolution of the measurements. Figure 10 shows the computed 
peak electric free field magnitudes for the single- and five-element source as a function of 
range along the center line. 

The extension of the single-dipole equations to a multiple source array by the method of 
superposition has proved to be a useful technique for prediction of the array environment. 

The computations have been supported by the measured data. The approach has demon- 
strated the amplitude degradation due to retarded time superposition; also, it has shown for 
the particular array dimensions under study in this section that the effect is beyond the sen- 
sitivity of the measurements. The discrepancies between the computed and measured data are 
therefore attributed mainly to the omission of mutual coupling losses from the computations. 
This argument is supported by the 27% higher computed value of II xT over the measured 
value of tlie component in figure 8. 

The theoretical and experimental approach has also shown that the free-field wave- 
shape radiated by the array can be obtained by measuring the tangential component of the 
magnetic field, H xT . This component interacts constructively with the ground and preserves 
the waveshape characteristics of the free field. The remaining principal components H and 
E , T interact destructively with the ground and therefore exhibit the sharp-spike waveform 
due to interference from the ground-reflected wave within a few nanoseconds after the wave- 
front arrival at the observation point. 
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Figure 10. Cofflouted electric free field for o single- ond five- 
dipole configuration. 


The results of the experimental and theoretical studies discussed in this section indicate 
that the array technique can indeed be used to increase Held strengths radiated by EMP simu- 
lators. The basic limitations of this technique have been identified to be: (1) distributed source 
retarded time effects and (2) mutual coupling losses. These limitations warrant further study 
and are discussed in the remaining sections of this report. . 


■ i i 

4. COLLINEAR DIPOLE ARRAY 

The collinear dipole array configuration shown in figures 11 and 12 was c onstructed to 
study the EMP environment of an array with a large aperture size compared w ope rise time 
of a free- field pulse. Two collinear banks containing three antenna^ ir each ferm the array. 

Each antenna has a current equivalent to the single-dipole current discussed in section 2. The 
z dimension is in the principal electric-field direction and corresponds to. more than ten times 
the free field pulse rise time of a single element. In this case, the different arrival times of 
the wavefronts from antenna banks A and B are clearly visible in figures 13 and 14, which 
show data obtained at observation points C and F In the test area. The initial peak in the 
waveforms is due to fields arriving from antenna bank A. The second peak is ^ue to the fields 
from antenna bank B which are displaced in time by a delay corresponding to the path length dif- 
ference between the observation point and points A and B. The path lengths are given in figure 11. 
A comparison of transit time differences agrees closely with the computed and measured results 
in figures 13 and 14. 

The field pattern of the collinear dipole array differs distinctly from the field patterns of 
the single dipole and the broadside array configurations. The polarity of the fields from each 
bank of dipoles in the collinear array is such that the nonprincipal field components (E^ , E_ t , 
H it ) tend to cancel and the principal field components (H^ , H , E f ) tend to add construe-' 
tively in the rectangular center region extending down range between A and B. Exact addition 
and cancellation of the principal and nonprincipal field components, respectively, occur only 
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Figura 11- Collinear dipole array geometry. 


on the centerline of the test area, defined along the line DG In figure 11. At other points in the 
rectangular center region, the largest contributions to the EMP environment from each anten- 
• a bank are due to the 9 independent H xT and H yT given by equations (12) and (13), and to the 
dominant sin 9 terms ,in equation (1*6). These terms have identical polarity for all antennas in 
the array and therefore add constructively in the field-equation summations. Because of this 
behavior, the collinear array has the desired characteristic of reinforcing principal field 
components and reducing nonprincipal field components within the center region. Ideally, this 
should increase the area of coverage by the principal fields compared with the coverage by the 
single dipole and the broadside array. The results, however, show that the gains in this re- 
spect are relatively small for the array dimensions under study. 
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Figure 12. The collineor dipole array. 
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The greatest limitation of the collineor array is the large time difference of wavefront 
arrival from the two antenna banks, which produce the two spikes in the waveforms. To 
eliminate this effect, the range of observation roust be increased to a point where transit- 
time differences are less than one rise time of the free-lield pulse from one antenna. Delta 
in figures 13 and 14 illustrate the decrease in transit-time differences between the 40- and 
the 180-ft ranges. Even at the far range, the waveshape degradation has not been eliminated. 
Points F and G would have to be at least 700 ft away from the array to absorb retarded time 
effects within a single-element pulse rise time at positions between F and G. 


5. EFFECTS OF MUTUAL COUPLING 

. i 

The previous sections predict the dipole-array environments under ideal lossless condi- 
tions. The approach was successful in providing early time waveshape characteristics and 
relative amplitudes of the EMP fields. One would expect that the lossless case provides an 
upper bound to the performance of a dipole array. It is of interest, however, to know the 
extent of losses caused by mutual coupling between array elements. Such information would 
lend more confidence in equations (12) through (16) for use in the design of larger array 
systems. The lossless case would give an upper bound to the environment, and the extent 
of mutual coupling would provide an efficiency factor for the design. 

! 

5.1 Measurements of Antenna Input impedance 

1 ' 

The network representation of mutual coupling between antennas Will be used in ^his study 
The approach is based on the knowledge of self and mutual impedance associated with the an- 
tennas. The antenna input impedance Z is measured by time-domain reflectometry (TDR). 
Figure 17 shows the measured input impedances for antenna elements I, ID, and V of the 
broadside array. In each case, the other four antennas weife removed for the measurements. 
The procedure used to make the measurement was to transmit a step-function input pulse 
through a 50- 0 reference cable to the antenna, and to observe the relative magnitude of the 
reflection (due to impedance mismatch) on a Tektronix 7704 oscilloscope. Standard TDR pro- 
cedure is used to compute first the reflection coefficient p at the junction of the reference 
cable and the dipole antenna, and then to obtain the antenna input impedance Z from: 

! , 

2 - R ohms (17) 

where R is the reference 50— 0 impedance and the value of p is giver, next to each oscilloscope 
picture in figures 17 and 18. 

Figure 17 shows the measured dipole antenna impedance as a function of antenna height 
above the ground. These measurements were taken to test the sensitivity of the TDR tech- 
nique. The results show that the input impedance Z increases a$ the antenna height above the 
ground is increased. This trend is correct if the antenna and its image in the ground are 
considered as a two-wire transmission line with characteristic impedance Z 0 given by 1 


^Romo, 5., Whirmpry, J. R., Fitldsond Wove® in Modern Rodio. Second Edition, John Wiley & Sone, Inc., N.Y. 
1964, p. 364. _ 
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Figure 17. Time-domain reflectometry. Dipole antenna impedance as a function 
of height above the ground. 


Z Q = 120 cosh -1 ( gj. 


(18) 


where d is the diameter of the two conductors and s is the separation distance. This equation 
applies for a line separated with an air dielectric and gives two times the impedance value for 
the antenna above a lossy ground plane. For a perfectly conducting ground plane, s is equal 
to twice the antenna height above the ground. The characteristic impedance Z 0 can then be 
computed from equation (18), assuming the perfectly conducting ground and an antenna diam- 
eter of 4 in. The results are compared with the measured data in table I. It is noted that the 
computed values of Z 0 are 1.5 to 1.8 times greater than the measured values of Z. The 
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Figure 18- Time-domain reflectometry. The effect of ad|acent elements on the 
measured impedance of the dipole antenna. 


trends are corrected, however, as evidenced by the comparison of relative impedance changes 
as a function of change in antenna height for both the computed and measured values given in 
table lb- The correlation of results is good, adding to confidence in the TDR measurements. 


5-2 Determination of Mutual Coupling Losses 

To study the effects of mutual coupling, the center antenna III is chosen as a typical ele- 
ment in the broadside array and the network representation is used to obtain the current and 
the electric field radiated by antenna in in the presence of adjacent elements II and IV. The 
mutual coupling losses introduced by the adjacent array elements can then be determined by 
comparison with performance of antenna in when all other antennas are removed. Measure- 
ments of impedance for antenna III were also taken with ail five elements present in the array. 
Mo difference in impedance was observed when elements I and V were removed; therefore, only 
the effect of adjacent elements on the performance of an antenna in the array is considered. 

Figure 18 shows the results of TDR measurements for the input impedance of antenna III 
in the presence of shorted and open circuited adjacent elements. In the open circuit condition, 
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Table la. Comparison of measured and commuted input impedance for a dipole 
antenna above a ground plane 


Antenna 

Number 

Antenna Height 
(Meters) 

Computed Z 0 
(ohms) 

TDR Measured Z 
(ohms) 

I 

3.96 

606 

366 

in 

2.44 

548 

320 

1 

0.91 

430 

284 


Table lb. Comparison between computed and measured relative impedance 
change as a function of change in antenna height 


Computed 

Measured 

Z(I)/Z(IH) = 606/548 = 1.11 
Z(m)/Z(V) = 548/430 = 1.27 
Z(I)/Z(V) = 606/430 = 1.41 

366/320 = 1.14 
320/284 = 1.13 
366/284 = 1.29 


the input impedance is the same as in the previous case where no other elements were present. 
This agrees with discussions of mutual coupling in antenna texts, namely that open-circuited 
elements should not effect the performance of active elements in an array. In the short-cir- 
cuit condition, the impedance changes from 320 to 254 0 indicating that the adjacent elements 
have a significant effect on the performance of antenna in. The equivalent circuit for this 
condition is shown in figure 19, where e, represents the voltage at the terminals of antenna 
in, and Z u and Z 2; represent the self-impedances of antenna in and the adjacent elements, 
respectively. The mutual impedance is represented by Z, 2 . The associated network relations 
are given by equations (19) and (20). 


= h z n-h z i2 < 19 > 

0 = - Ij Z 12 + I 2 Z 22 . (20) 


Active 

Antenna 


Z n _Z 12 Z 22* Z 12 



Passive 

Shorted 

Adjacent 

Elements 
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The mutual impedance is determined experimentally by the following procedures: 3 

(1) The self-impedance Z n of antenna in is measured with all other elements removed. 

(2) The impedance Z 1 into antenna in is measured with shorted adjacent elements 3 
and IV. 

The mutual impedance is then computed from; 


^12 = v z u ( Z ll ~ Z jl' 


( 21 ) 


Equation (21) assumes adjacent shorted elements with Z u = Z for identical antenna ele- 
ments. One obtains Z u and Z t from figures 17 and 18, respectively, and the mutual impedance 
Z 12 can be computed from equation (21): 


Z 12 = /320 (320 - 254) = 146 H 


(22) 


Next, determine the voltage e , across the biconic section of antenna III as follows. The gen- 
erator voltage V ( is 30 kV at the pulser. This in turn is transmitted through a 100-fi line to 
the input at the biconic section of the antenna. At this point there is a mismatch in impedance 
with an associated reflection coefficient f : 


P _ Z - 100 
Z + 100 


(2a 


where Z = 320 fi, the self-impedance of antenna in. Compute e. by: 


e, = (l-r)v f (24) 

In this case, r = 0.52 and e, = 14.4 kV. Sufficient data are available to compute the currents 
in equations (19) and (20). The computation results in I, = 56.83 amperes (A), and I 2 = -25.9 A. 
Therefore, in the short-circuited adjacent element condition, antenna in induces a net current 
of -25.9 A in the adjacent elements. This current in turn radiates an opposing field to the field 
of the active element III and causes a decrease in resultant field strength. 

A formula for computing the peak free space electric-field strength at r on the center- 
line of a bicone-fed dipole antenna is: 3 


60 v _ 60 I 

TT ’ ~ 


%1anten, R. C. Microwave Scanning Anfanna*, Volume II Array Theory ond Practice. Academic Prew, N.Y. 
op. 157-160. 

3 Schelkunoff, S. A., Friii, H. T., Antennal, Theory ond Practice. John Wiley ond Soni, Inc., N.Y. 1952, p. 
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(25) 

1966, 

106. 




where v is the terminal voltage and Z is the self-impedance of the antenna. In our case, 
v = e 1 = 14.4 kV, Z = 320 fi , and r = 12.2m. The net antenna current I is equal to 45A in the 
absence of other antenna elements. This corresponds to a free- field strength of E Q = 222 
V/m for antenna ID. In the lossy case, where the adjacent short-circuited elements reduce 
the net currunt to I = 56 <83 - 25.9 = 30.93 A the corresponding field strength is E 0 = 152 V/m. 
(In the strict sense, the subtraction should include retarded time effects. However, the ele- 
ment spacing is so close that a direct subtraction is justified.) Consequently, mutual coupling 
has reduced the efficiency of antenna III to 69% of its original output. 

The mutual coupling studies have shown that significant losses can arise from parasitic 
adjacent elements in an array. Therefore, design considerations should include these effects. 
The arguments presented in this section have been based on the impedance measurements and 
the network representations of mutual coupling. From these considerations, we obtained the 
net current and field strengths required to obtain the decrease in efficiency of a radiating ele- 
ment in the presence of adjacent antennas. In the following section, the efficiency of the arrays 
studied in this report is discussed from the standpoint of measured fields. Thus, one has the 
results of two approaches to estimating array efficiency: one approach dependent upon knowl- 
edge of Impedances and voltages in the array (and independent of the radiated fields), and the 
other approach dependent upon knowledge of the radiated fields (and independent of the imped- 
ances and voltages in the array). 


6. EFFICIENCY FACTOR OF THE ARRAY MODELS 


Two mechanisms for degradation of performance have been identified which introdr -e 
limitations in the dipole technique for EMP simulation. Neglecting other losses, the physical 
size of the array determines a minimum range where the waveshape is not degraded by the 
arrival time differences from extreme elements. The array is a distributed source for prac- 
tical ranges; therefore, the field strengths do not reach the ideal N times the single-element 
output level, where N is the number of antennas in the array. 

Define a lossless efficiency factor (E.F.) that provides a measure of bow much the dis- 
tributed nature of the array affects the output, cs compared with the output of a point source 
having N times the single-array element strength. The H xT component is used as a reference 
field for the computations. 


LOSSLESS E. F. = 


'computed H xT for lossless array 
^computed H xT for single element 


x 100% 


(26) 


The E.F. factor approaches 100% at ranges where the array approximates a point source. 
Table II shows computed values of the lossless E.F. for the array configurations studied in 
this report. The efficiency computed from H is, in general, less than that obtained from 
the free-field output, as evidenced by the results shown for the five-element array in figure 
10. For that case, the free-field efficiency approached 100%, whereas the use of H xT in equa- 
tion (26) has resulted in an efficiency factor of 76%. 

Additional losses are introduced by mutual coupling effects. These losses are included in 
a second E.F. which iB referenced to the measured H xT component, and which includes any 
other real losses present in the array. 
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ANTENNA 

CONFIGURATION 


Table n. The array efficiency factor 
REFERENCE FIELDS (A/m) 


LOSSLESS 
H T 


MEASURED 
H „ 


EFFICIENCY FACTOR 


LOSSLESS 

MEASURED 

E.F. 

E.F. 

1001 

67* 

100% 

88* 

99* 

84* 

76* 

59* 


/ i\ / measured array H_. \ 

MEASURED E. F. =(I] II * 100% 

\N/\measured singl e element H _/ 


Results of the computations are given in table n. 

For antenna m with short-circuited adjacent elements, the measured E.F. is 67%, This 
agrees closely with the 69% efficiency predicted by the mutual coupling results. In the open 

















I 



circuit adjacent element condition, the measured E.F. is 88%. This indicates that losses are 
introduced which were not apparent in the impedance measurements. The early time fields 
from the collinear array can be considered as a three-element broadside array. The meas- 
ured E.F. is 84% for this condition. By comparison, the five-element broadside array measured 
efficiency is 59%. This indicates that increasing the number of elements in order to increase 
field strength has a practical limit, beyond which little is gained by adding elements. 


7. CONCLUSIONS 

, i 

This report summarizes the results of studies conducted at the HDL Electromagnetic Ef- 
fects Laboratory on nuclear electromagnetic pulse simulation by dipole antenna array tech- 
niques. It i3 evident from these studies that the aperture size and range of operation are im- 
portant factors that must be taken into consideration in the dipole array design. The dimen- 
sions of interest are the path length difference between extreme elements in the array and the 
point of observation in the test area. The transit-time differences must be at least smaller 
than the rise time of the single-element free space pulse to prevent vavesha!pe degradation. 
Large apertures in either or both the electric and magnetic directions may cause this . 
degradation. 

Losses from mutual coupling must also be considered In design procedures. The method 
presented in this report can be used to estimate these losses by taking impedance measure-' 
mants on reduced scale models, or possibly on full size models of large antenna elements. 

The efficiencies computed for the array configurations studied in this report should be con- 
sidered as representative examples of the correlation between the various impedance, current, 
voltage, and field parameters for dipole antenna arrays. The efficiency factors indicate that 
optimal arrays should have as few antenna elements as possible to reduce losses introduced by 
mutual coupling and by distributed source effects. The theoretical methods and the computer 
programs developed during this study have been supported by experimental results, and can 
be used to predict the output of such arrm™ 
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